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摘要 :水 资源 是 保障 我 国 西北 干旱 半 干 旱地 区 生态 环境 安全 的 关键 因素 。 以 新 疆 奎屯 河流 域 为 例 , 通 过 修正 SCS 模型 土壤 持 水 
量 及 初 损 率 参数 计算 方法 ,寻找 适用 于 干旱 半 干 旱地 区 山区 典型 流域 春季 融雪 期 径流 模拟 模型 ,为 流域 掌握 水 资源 量 及 生态 用 
水 提供 决策 依据 。 与 以 往 研究 不 同 之 处 在 于 :首先 ,引入 度 -日 模型 修正 降水 量 参数 ,以 满足 流域 降雨 -融雪 混合 补给 径流 特征 。 
其 次 ,利用 多 期 MODIS 数据 驱动 的 Ts/VI 特征 空间 理论 结合 土壤 水 分 吸收 平衡 原理 计算 土壤 持 水 量 参 数 (S) ;再 运用 聚 类 分 析 
法 对 初 损 率 (A) 取 值 方法 进行 改进 。 通 过 参数 算法 改进 后 的 SCS 模型 ,参数 率 定期 和 验证 期 纳什 效率 系数 和 相对 误差 系数 分 
别 为 0.92 和 0.64,0.7% 和 -1.5% 。 结 果 表 明 :1) 参 数 算法 改进 后 SCS 模型 能 实现 奎屯 河流 域 春季 融雪 期 日 径流 模拟 。2) 利用 
遥感 大 尺度 地 表 信息 参数 化 技术 反 演 SCS 模型 参数 ,实现 了 遥感 数据 为 SCS 模型 提供 大 尺度 空间 数据 的 同时 ,间接 实现 了 模 
型 参数 由 点 状 数据 向 面 状 数据 转化 的 可 能 ;3) 初 损 率 (入 ) 多 组 取 值 法 可 有 效 提高 干旱 半 干 旱地 区 大 尺度 流域 径流 模拟 精度 。 
关键 词 :奎屯 河流 域 ;融雪 期 ;SCS-CN 模型 ;参数 算法 改进 ;遥感 
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Abstract: Water resources)are key factors of ecological environmental security in northwest arid region of China. They are 
also the, mostMimportant factors for socio-economic development against the background of global warmer, especially in arid 
regions. Tis necessary for arid regions to calculate total water resources because it can provide a reference for the 
governnient with which to formulate strategies. Water resources may have a large area and be supplied by runoff from 
mountain snowmelt and precipitation. The goal of this paper was to determine the suitable method to simulate runoff in arid 
areas., The Soil Conservation Service Curve Number (SCS-CN) developed by the U.S. Department of Agriculture National 
Resources Conversion Service ( NRCS) is the most popular and widely applied model for direct runoff estimation. This 


method was modified by accounting for the static portion of infiltration and the antecedent moisture. This model has 
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stimulated a great deal of discussion among scientists and hydrologists. The model is based on the water balance equation 
and curve number CN, which is derived from the tables given in the National Engineering Handbook for catchment 
characteristics such as soil texture, land use, hydrologic condition, and initial soil moisture condition. Based on the 
spatiotemporal differences among watersheds, international and domestic academics have developed different methods to 
improve the SCS-CN model. One option is to improve its mechanism and another is to improve the parameter calculation 
methods. Because there is considerable scope to improve the SCS-CN model, we discuss a parameter algorithm to improve 
the method for snowmelt and precipitation mix and large-scale basins in arid regions as a solution to a major problema THis 
study focused on the Kuitun River Valley. We explored the adoption of an SCS model runoff simulation in arid and semi-arid 
regions with snow-melt and rainfall in spring by modifying the calculation method of SCS model parameters. To satisfy the 
characteristic of mix supplied runoff, precipitation was revised to represent the sum of rainfall and snowmeltt The snowmelt 
was calculated by the degree-day model. This was the first time MODIS satellite products with approkimately Wm resolution 
were used to invert the Land Surface Temperature and Normalized Difference Vegetation Index.SLhen¥ we used the surface 
temperature/ vegetation index (Ts/VI) constructed in a 2D scatter plot. The combined soil moisture absorption balance 
principle was used to calculate the moisture-holding capacity of soil. We used clusten analysis to modify the initial 
abstraction computing methods. The calibration and validation periods of Nash-Sutcliffe Mefficieney were 0.92 and 0.64, 
respectively. Relative errors were 0.7% and -1.3%, respectively. This indicated thatthe/ improved model was effective in 
simulating spring runoff in the Kuitun River Valley. Using remote sensingsparameter information technology to improve the 
SCS model can indirectly implement data conversion from point to plane. WEstablishing a database of the initial abstraction 
can improve the precision in effectively simulating runoff in large-s€ale basins in arid and semi-arid regions. To circumvent 
the bottleneck caused by lack of data, reference to simulated ruiioff ean be used under similar basin conditions in data- 


lacking regions. 


Key Words: Kuitun River basin; snowmelt period SES:CN model; improved algorithm of parameters ; remote sensing 


中 国 西北 部 干旱 半 干 旱地 区 河流 年 要 发 源 于 山区 ,地 表 径 流 主要 由 高 山 带 冰 ( 川 ) 雪 融 水 、 中 山 森 林带 降 
水 和 低 山 带 基 岩 裂 隙 水 等 组 成 。 长 期 依靠 自然 界 独 特 的 水 循环 过 程 维 持 着 脆弱 的 平衡 关系 "1 。 在 全 球 气候 
变 暧 影响 下 ,干旱 区 内 陆 河 流域 水 资源 的 不 确定 性 不 断 加 剧 ,直接 影响 流域 生态 安全 。 尤 其 在 山区 及 其 周边 
流域 ,每 年 4 一 5 月 融雪 期 $ 经 流 时 空 分 布 差异 显著 ,频繁 发 生 的 融雪 型 洪水 ,使 流域 大 面积 植被 遭 到 破坏 ,水 
土 流失 严重 ,生态 环境 十 分 脆弱 ,给 当地 农业 生产 和 人 民生 活 带 来 巨大 安全 隐患 。 因 此 ,定量 研究 山区 流域 地 
表 径 流 ,对 认识 流域 水 资源 季节 性 变化 规律 ,开展 生态 环境 保护 和 防洪 减灾 工作 具有 重要 意义 。 

在 流域 地 表 径 流 模拟 过 程 中 ,SCS-CN 模型 ( Soil Conservation Service Curve Number ,简称 SCS 模型 ) ,是 目 
前 全 球 应 用 最 为 六 泛 的 降雨 -径流 模型 之 一 。 由 原 美国 农业 部 水 土 保持 局 (USDA-SCS ) 根据 境内 不 同 地 区 流 
域 降雨 -径流 资料 的 搜集 、 整 理 .分析 人 研究 得 出 的 经 验 模型 。 具 有 结构 简单 物理 概念 明确 、 所 需 参 数 少 等 特 
点 , 因 被 SWAT、Mike Hyro Basin 等 分 布 式 水 分 模型 应 用 于 地 表 径 流 模拟 而 在 世界 范围 内 得 到 广泛 应 用 。 该 
模型 宝 要 研究 集中 于 湿润 地 区 ,上 且 深入 研究 者 甚 多 “5 。 王 旱 半 干旱 地 区 也 有 应 用 ,例如 周 淑 梅 .张秀英 及 李 
舟 等 人 通过 对 模型 的 改进 ,分 别 在 陕西 .甘肃 和 西北 高 寒山 区 小 流域 模拟 都 取得 了 良好 的 效果 “1 。 地 域 广 
阔 的 干旱 区 ,流域 尺度 大 ,具有 降雨 冰雪 融 水 混合 补给 径流 等 特征 ,其 地 表 径 流 模 拟 区 别 于 单纯 以 降雨 驱动 
的 湿润 地 区 小 流域 径流 模拟 ;加 之 山区 可 达 性 差 , 人 为 观测 极为 困难 ,利用 遥感 大 斥 度 信息 观测 优势 结合 水 文 
模型 ,寻找 适合 于 干旱 半 干 旱 区 地 表 径 流 模拟 方法 ,对 当地 生态 环境 、 农 业 生 产 及 人 类 生活 具有 重要 现实 意 
义 。 基 于 此 ,本文 预 借助 遥感 技术 对 SCS 模型 参数 计算 方法 进行 改进 以 提高 径流 模拟 精度 ,具体 改进 方法 如 
下 :1) 借 助 遥 感 大 尺度 地 表 空 间 信 息 参 数 化 特性 ,结合 MODIS 卫星 数据 驱动 的 特征 空间 理论 反 演 流域 空间 土 
壤 含 水 量 , 再 利用 土壤 水 分 吸收 平衡 原理 ,间接 实现 模型 土壤 持 水 量 ($) 空 间 数据 的 估算 ;2) 针对 流域 降雨 、 


http://www.ecologica.cn 


201707.00097V1 


图 
国 


XIV 


hina 


ChinaXiv 合 作 期 刊 


4458 生 态 学 报 37 卷 


融雪 混合 补给 径流 特征 ,引入 度 - 日 模型 ,将 降水 量 参数 修正 为 降雨 与 积 雪 消 融 水 当量 之 和 ;3) 采 用 聚 类 分 析 
法 将 参数 率 定 期 初 损 率 (A ) 进行 聚 类 ,选择 多 组 值 进行 模拟 ,以 提高 模拟 精度 。 小 流域 SCS 模型 应 用 时 汇 
流 时 间 较 短 可 忽略 汇流 项 ,大 尺度 流域 SCS 模型 应 用 要 确保 汇流 时 间 小 于 24 h。 


1 数据 与 方法 


1.1 研究 区 概况 

本 文 以 新 疆 天 山 文 脉 依 连 哈 比尔 泰山 北 划 奎屯 河 
流域 为 研究 区 (图 1) ,流域 面积 1945km ,位 于 83。 
30' 一 85°08”E,43°30' 一 45°00”N 之 间 , 是 天 山北 坡 中 
段 仅 次 于 玛 纳 斯 河流 域 的 第 二 大 流域 。 集 水 区 面积 较 
大 ,年 平均 气温 9%C ,海拔 为 1121 一 4909m , 其 中 海拔 
3700m 以 上 为 冰川 及 永久 性 积 雪 有 覆盖 ,致使 流域 径流 补 
给 多 样 , 既 有 高 山 冰 川 和 永久 性 积 雪 融 水 补给 ,又 有 中 


44°20'00"N 


44°00'00” 


低 山 区 降雨 和 季节 性 冰雪 融 水 补给 。 每 年 4 一 5 月 为 流 § 和 : 
域 春季 典型 融雪 期 ,径流 补给 以 积 雪 融 水 为 主 ,降雨 为 有 < 
辅 ,, 当 气温 又 然 上 升 时 会 引发 季节 性 融雪 洪水 。 基 于 So ey 0 
以 上 特征 ,本 文选 取 流 域 出 山口 将 军 庙 水 文 站 2005 一 图 1 奎屯 河流 域 示意 图 

2007 年 日 观测 数据 ,模拟 奎屯 河 出 山口 4 一 5 月 地 表 径 Figl Watershed location of Kuitun river basin 
流量 ,为 周边 灌区 及 城市 的 防洪 减灾 和 水 资源 分 配 提供 

参考 依据 。 


1.2 数据 资料 

奎屯 河流 域 水 文 模拟 需要 的 基础 地 理 数据 包括 :气象 水文 .卫星 遥感 影像 等 多 源 数据 类 型 。 其 描述 和 来 
源 见 表 1。 其 中 ,为 获取 人 研究 区 参数 率 定期 地 表 温 度 及 植被 指数 月 均值 数据 , 需 下 载 空 间 分 辨 率 为 1 km 的 
MODIS 陆地 3 级 标准 数据 产品 MOD13A3W MOD11A2 , 并 对 数据 进行 预 处 理 , 包 括 对 数据 进行 投影 转换 、 几 
何 校正 和 矢量 边界 裁剪 等 ,得 到 覆盖 全 区 的 栅 格 数据 。 已 有 大 量 人 研究 证 实 MODIS 卫星 产品 精度 能 有 效 反 映 
区 域 地 表 温 度 及 植被 指数 情况 0 。 


表 1 奎屯 河流 域 基础 地 理 数 据 


Table 1 Geographical data for kuitun river basin 


数据 类 型 Data type 数据 描述 Data description 数据 来 源 Data sources 
新 疆 维 吾 尔 治 区 局 观测 
气象 Meteorological 奎屯 河流 域 2005 一 2007 年 4 一 5 月 逐日 降水 气温 观测 数据 维吾尔 明治 区 气象 局 观测 
贷 锌 
水 交 志 中 量 奎屯 河流 域 出 山口 将 军 庙 水 文 站 2005 一 2007 年 4 一 5 月 日 径流 ”新 疆 维吾尔 自治 区 水 文 资源 局 观 
数据 测 资料 
i 本 野外 实测 将 军 庙 水 文 站 周边 13 个 
奎 市 河 ! 证 质 Se 下 用 外 实 泊 人 二 
土壤 求 分 Soil moisture ee 2006 年 4 一 5 月 野外 实测 土壤 全 0 一 10cm 土壤 剖面 样本 采集 ,实验 
室 分 析 
2005 一 2006 年 4 月 .5 月 归 一 化 植被 指数 (NDVI) 数据 来 源 于 
MODIS 陆地 3 级 标准 数据 产品 MOD13A3 ,空间 分 辩 率 lkm, 月 尺度 
卫星 遥感 影像 数据 MODIS 产品 来 源 于 NASA 网 站 
Remote sensing image 2005 一 2006 年 4 月 .5 月 地 表 温 度 (LST) 数据 来 源 于 MODIS 陆地 3 DEM 数据 来 源 于 地 理 空 间 数 据 云 


级 标准 数据 产品 MOD11A2 ,空间 分 辨 率 lkm, 月 尺度 数据 奎屯 河流 
域 SRTM30m 分 辨 率 DEM 数据 
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2 研究 方法 


2.1 SCS 模型 原理 
SCS 模型 是 基于 水 量 平衡 方程 (1) 存 在 的 : 
P=lat+F+0 (1) 
式 中 ,P 为 总 降雨 量 (mm) ;1 为 地 表 径流 生成 前 降雨 量 初 损 值 ,包括 地 面 填 洼 截留 表层 鞭 水 和 下 渗 的 初 损 
量 (mm) ;F 为 地 表 径流 生成 后 的 后 损 (mm) , 即 实际 累积 下 渗 量 (不 包括 7,) ;0 为 地 表 径 流量 (mm) 。 
与 此 同时 ,Mockus 等 人 通过 大 量 实测 数据 分 析 建 立 了 降雨 -径流 关系 式 : 


0 -地 
P-la 5S (2) 


式 中 ,5 为 土壤 当时 最 大 可 能 持 水 量 ,是 后 损 的 上 限 (mm)。 由 于 式 中 7 数据 不 易 获取 ,通常 引入 参数 初 损 率 
(和 A) ,建立 1 与 5 的 线性 关系 式 : 


Ja=AS (3) 
式 中 ,A 为 初 损 率 ,具有 区 域 化 特征 。 将 方程 (2) (3) 带 入 (1) 中 得 SCS 径流 模拟 核 必 方程 为 : 
0 (P>A5, 盏 则 0=0) (4) 


式 中 ,0 为 地 表 径流 深 (mm) ;P 为 总 降水 量 (mm) ;5 为 当时 最 大 可 能 持 水 量 (mm) 。A 为 初 损 率 。 其 中 $ 值 
通常 由 CN 计算 而 得 ,计算 方法 如 方程 (5) 所 示 : 
Se= 25400 a 国际 制 单位 ) (5) 

式 中 ,CN 值 为 径流 曲线 数 ,通常 由 土地 利用 方式 .土壤 质地 和 和 降雨 事件 前 期 土壤 湿润 情况 (antecedent moisture 
condition, AMC) 等 数据 来 确定 。 其 详细 分 类 标准 圳 参考 美国 (国家 工程 手册 》'"。 

初 损 率 (A) 通 常 取 标 准 值 0.2 进行 计算 ;但 直 于 不 同 流域 具有 一 定 的 时 空 差异 ,A 值 取 标 准 值 不 能 满足 不 
同 流域 的 需求 ,本 文 利用 实测 降雨 .径流 数据 和 计算 s 值 结合 公式 (6) 对 初 损 率 ( 和) 进行 参数 率 定 ,具体 方法 
如 下 : 


me 2P, - 0, 一 /0 ° + 401S 
2S 

式 中 ,Pi、0, 分 别 为 模型 参数 率 定期 流域 实测 日 降水 量 (mm) 和 日 径流 量 (mm) 。 
2.2 SCS 模型 参数 计算 方法 的 改进 
2.2.1 降水 量 计算 方法 的 改进 

由 于 原 有 SCS 模型 对 降雨 量具 有 一 定 条 件 限 制 , 即 模型 认为 当 降 雨量 (P) 大 于 AS 时 ,产生 地 表 径 流 ,和 否 
则 地 表 径 流量 为 零 。 该 条 件 使 SCS 模型 无 法 应 用 于 降雨 量 稀少 的 干旱 半 干 旱地 区 。 本 文 根 据 研究 区 4 一 5 月 
降雨 和 积 委 融 水 混合 补给 径流 的 特点 ,将 SCS 模型 降水 量 修正 为 流域 降雨 量 与 积 雪 消 融 水 当量 之 和 , 如 方程 
(7) 鼎 示 。 


(6) 


P=P,+P. (7) 
式 中 ,P 为 SCS 模型 径流 方程 输入 的 总 降水 量 (mm)。P, 为 实测 日 降雨 量 (mm);P, 为 日 积 雪 消融 水 当量 
(mm) ,引入 度 日 模型 进行 计算 。 
度 -日 模型 是 基于 冰雪 表面 温度 建立 的 冰雪 融 水 当量 计算 模型 。 广 泛 应 用 于 北欧 .阿尔卑斯 山 青藏 高 原 
等 地 区 的 冰雪 消融 研究 中 。 优 点 在 于 模型 计算 简单 参数 容易 获取 且 计 算 精 度 较 高 。 因 此 ,本 文 引 入 度 -日 模 
型 计算 研究 区 4 一 5 月 积 雪 消 融 水 当量 ,具体 计算 方法 如 下 : 
P, =DDF x PDD (8) 
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式 中 ,P 为 日 积 雪 消 融 的 水 当量 (mm) ;DDF 为 冰川 或 雪 的 度 日 因子 (mm qd"'%C "1) ,利用 张 勇 5531 等 人 根据 海拔 
对 中 国 西北 地 区 冰雪 度 日 因子 计算 结果 ,采用 克 里 格 空间 插值 法 计算 而 得 ;PDD 为 某 一 时 段 内 的 正 积 温 ， 
PDD 可 通过 下 式 计算 : 


PDD= SHxT, (9) 
t=1 


式 中 ,7 为 某 时 段 的 平均 气温 (%C ) ; 太 , 为 逻辑 变量 , 当 7,=0%C 时 ,有 H=1.0; 当 7 <0C 时 ,有 H,=0.0。 
2.2.2 土壤 持 水 量 (5) 计 算 方 法 的 改进 

原 有 SCS 模型 $ 值 是 通过 CN 值 进行 计算 的 ,存在 涉及 参数 多 ,数据 难以 获取 有 旦 计算 过 程 复杂 等 问题 。 
与 传统 方法 相 比 ,本 文 改进 之 处 在 于 ,首次 利用 遥感 及 7.-VI 特征 空间 理论 反 演 土 壤 含 水 量 参 数 结合 土壤 水 
分 吸收 平衡 原理 计算 $ 值 ' 中 ,计算 方法 如 下 :， 

S= Ws —W.o (10) 

式 中 ,5 值 为 当时 最 大 可 能 持 水 量 ; WW 为 土壤 饱和 含水 量 , 通 过 测定 野外 土壤 样品 的 汪 二 隆 际 度 和 土壤 容重 
计算 而 得 ,由 于 数据 缺乏 ,本文 取 最 大 值 100%'' ;WW 为 平均 土壤 含水 量 ; 本文 下 ,为 研究 区 参数 率 定期 
2005 一 2006 年 4 一 5 月 土壤 含水 量 月 平均 值 ,是 利用 MODIS 卫星 数据 产品 地 表 温 度 人 7 ) 和 归 一 化 植被 指数 
(NDVI) 构 建 的 温度 -植被 干旱 指数 (7T-VI) 特征 空间 理论 进行 估算 的 上 为 使 SCS 模型 中 参数 单位 统一 为 mm， 
需 用 水 层 厚度 将 土壤 含水 率 转换 为 土壤 含水 量 ,mm ,进行 计算 。 

Moran' 等 人 在 地 表 温 度 ( 7 ) .植被 指数 (WI) 与 土壤 湿度 (Water Gontent of soil, SRWC ) 之 间 关 系 研究 中 
发 现 , 土 壤 湿 度 与 二 者 变化 具有 极为 密切 的 关系 , 若 将 区 域内 每 个 像 元 的 下 值 作为 x 轴 ,7 作为 y 轴 ,可 构建 
出 呈现 阶梯 状 的 特征 空间 二 维 散 点 图 , 散 点 图 的 两 条 边界 线 分 别 代表 区 域 土壤 湿润 程度 的 干 边 和 湿 边 ,再 通 
过 干 , 湿 边 方程 系数 结合 NDVI 计算 得 到 温度 植被 干旱 指数 TVDI 值 , 该 值 可 直接 反映 地 表土 壤 水 分 情况 , 计 
算得 到 任 一 点 的 TVDI 值 介 于 0 和 1 之 间 。TVDI 值 武大 ;代表 土壤 含水 量 越 低 ,反之 土壤 含水 量 则 越 高 。 

为 了 使 计算 结果 更 接近 实测 值 ,将 TVDK 值 与 实测 土壤 水 分 值 进 行 拟 合 修正 ,最 终 得 出 区 域 土壤 水 分 
(SRWC) 值 。 其 中 ,TVDI 计算 公式 如 下 : 


3 


TVDI= 11 
7 (11) 


式 中 ,7 和 了 ,分 别 表 示 区 域 相同 植被 指数 条 件 下 所 对 应 的 最 高 温度 和 最 低温 度 ,计算 公式 如 下 : 
T=aVi+t6b; 7 =cYI+di (12) 
式 中 ,ob cd 分别 为 干 齐 边 拟 合 方程 系数 。 
最 后 ,根据 公式 (13) Qu4) 求 算出 研究 区 土壤 含水 量 (SRWC)09 


SRWC=SRWCw -TVDI( SRWC, -SRWC,) (13) 
式 中 ,SRWG@, 为 湿 边 上 对 应 的 最 大 值 ,取水 体 湿 度 100%。SRWC, 为 干 边 上 对 应 的 最 小 值 ,计算 公式 为 : 
SRWC, = (100-(100-Y.)/X,)/100 (14) 


式 中 ,六 为 对 应 该 点 植被 全 和 覆盖 像 元 TVDI 值 ,7 为 对 应 该 点 的 实测 SRWC 值 。 
2.2.30 初 损 率 (A) 计算 方法 的 改进 

目前 ,已 有 的 和 值 算法 有 3 种 方法 。 即 标准 值 法 (A=0.2) .事件 分 析 法 和 反 算 法 , 反 算 法 应 用 较 多 , 即 : 每 
一 组 降水 ,径流 数据 结合 5S 值 可 计算 出 与 之 一 一 对 应 的 A 值 ,再 通过 计算 A 的 中 值 、. 众 数 .平均 值 等 数学 方法 
选取 模拟 精度 最 高 的 ,作为 最 终 唯一 的 和 A 取 值 。 

本 文 改 进 之 处 在 于 ,利用 聚 类 分 析 法 将 降水 ,径流 事件 结合 公式 (6) 计 算 所 得 的 若干 A 值 进行 聚 类 分 析 ， 
得 到 30 组 不 同 的 聚 类 中 心 值 ,并 将 这 30 组 聚 类 中 心 值 作为 流域 的 A 取 值 库 , 再 根据 迭代 次 数 计算 的 初始 中 
心 最 小 距离 将 所 有 和 值 及 其 对 应 的 降水 ,径流 事件 分 别 归 类 到 这 30 组 数据 中 ,使 每 组 降水 、 径 流 数据 都 有 与 
之 对 应 的 聚 类 中 心 A 值 ,该 方法 与 传统 和 取 唯 一 值 相 比 ,模拟 结果 更 接近 实际 值 。 
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3 ”参数 算法 改进 后 模型 参数 率 定 


3.1 降水 量 参 数 的 确定 

奎屯 河流 域 降水 量 为 降 南 量 和 积 雪 消 融 水 当量 之 和 ,结合 流域 2005 一 2007 年 4.5 月 实测 降雨 量 和 温度 
数据 ,计算 出 该 时 段 的 修正 降水 量 见 图 2, 图 中 绘制 了 183d 的 降水 量 修正 值 ,其 中 降水 量 为 0 值 代表 该 日 气 
温 在 0%C 以 下 且 无 降雨 。 此 外 ,由 于 积 雪 融 水 当量 的 主要 决定 因素 是 温度 ,而 研究 流域 的 日 平均 气温 从 省 月 
到 5 月 是 呈现 逐步 上 升 趋势 的 ,因此 所 有 年 份 修正 降水 量 的 分 布 都 表现 为 4 月 小 于 5 月 。 
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图 2 奎屯 河流 域 2005 一 2007 年 4\S 月 降水 量 修正 值 
Fig.2 Modified precipitation of Kuitun River Basin in April and May from 2005 to 2007 


3.2 S 值 的 确定 

采用 2005 一 2006 年 4.5 月 122d 日 数据 作为 SCS 模型 的 参数 率 定期 计算 当时 最 大 可 能 持 水 量 (S) 值 。 
3.2.1 流域 边界 的 确定 

利用 研究 区 30m 分 辨 率 DEM 人 和 数据 ,借助 Mike Hyro 
Basin 软件 ,生成 流域 边界 (图 8) 。 
3.2.2 ”地 表 温 度 及 植被 指数 计算 

按 流域 边界 提取 奎 电 河流 域 2005 一 2006 年 4 月 一 
5 月 共 32 期 地 表 温 度 和 归 一 化 植被 指数 ,生成 流域 
2005 一 2006 年 4 一 5 月 平均 地 表 温 度 空 间 分 布 图 和 归 


44°00'00”"N 


43°50'00” 


一 化 植被 指数 空间 分 布 图 (图 4, 图 5)。 入 
3.2.3 /pI 话 生 3 
) 计算 二 湿 边 方程 84°15'00”" 84°30'00”" 84°45'00"E 
oa 、 图 例 一 河流 闭合 流域 
利用 到 ,- 太 特征 空间 理论 ,结合 生成 的 2005 一 2006 Wm > CA 


年 MM 一 5 月 平均 地 表 温 度 和 平均 归 一 化 植被 指数 空间 分 图 3 MIKE Hyro Basin 生成 奎屯 河流 域 边界 
布 图 生成 二 维 散 点 图 和 拟 合 方程 ,并 结合 公式 (11) 计 Fig.3 Kuitun River Basin boundary generated by MIKE 
算出 研究 区 2005 一 2006 年 4 一 5 月 平均 TVDI 空间 分 布 ”YDRO BASIN 
图 (图 6) 。 

(2) 土壤 含水 量 计算 

由 于 流域 可 达 性 较 差 ,使 野外 采样 点 布设 受 限 , 现 有 野外 实测 数据 为 2005 .2006 年 4.5 月 研究 流域 将 军 庙 水 
文 站 周边 13 个 土壤 剖面 0 一 10cm 实测 土壤 含水 量 (SRWC) 值 , 样 点 布设 为 每 隔 1km 进行 一 次 采样 ,每 个 采样 点 
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采取 3 次 重复 采样 的 方法 ,将 重复 采样 土壤 样本 研磨 均匀 混合 装 铝 盒 后 称 重 , 带 回 实验 室 采 用 人 烘 干 法 测定 样本 
土壤 含水 率 。 然 后 根据 方程 (13) 计 算出 流域 2005 一 2006 年 4 一 5 月 平均 土壤 含水 量 空间 分 布 图 (图 7) 。 


尼 N 
S 入 
上 
mn 
be 
己 
呈 
ww 
HT 
守 | 图 例 
mm 高 : 301.38 
低 : 199.66 0 10 20 m 
二 1 1 
84°15'00” 84°30'00”" 84°45'00"E 


图 4 奎屯 河流 域 2005 一 2006 年 4 一 5 月 平均 地 表 温 度 空 间 分 布 
Fig.4 Spatial distribution of average land surface temperature in 


April and May from 2005 to 2006 in kuitun river basin 
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6 ”奎屯 河流 域 2005 一 2006 年 4 一 5 月 平均 TVDI 值 空间 分 布 
Fig.6 Spatial distribution of TVBI in April and May from 2005 


to 2006 in Kuitun River Basin 


关于 遥感 反 演 土壤 含水 量 精度 是 否 可 靠 , 采 用 计算 
值 与 实测 值 对 比 的 方法 进行 验证 。 即 用 时 外 实测 采样 
点 坐标 提 职 SRWQ 空间 分 布 图 中 对 应 的 值 ,并 与 实测 
土壤 含水 量 进 簿 相关 分 析 ( 图 8) ,得 相关 系数 + 为 0. 
7406, 相 对 误 荆 系数 RE 为 -1.3% ,该 结果 表明 ,遥感 反 
演 填 过 合 水 量 精度 可 靠 。 

利用 计算 所 得 流域 平均 土壤 含水 量 值 结合 公式 
(10) ,初步 确定 奎屯 河流 域 4 一 5 月 SCS 模型 参数 当时 
最 大 可 能 持 水 量 (S) 值 为 80.5mm。 
3.3” 初 损 率 (和) 的 确定 

根据 计算 所 得 $ 值 和 流域 参数 率 定期 及 验证 期 
183d 修正 降水 量 和 观测 径流 量 数据 结合 公式 (6) 计 算 
出 183 组 初 损 率 (A ) 值 ,以 模型 次 降水 量 大 于 下 渗 量 
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图 5 奎屯 河流 域 2005 一 2006 年 4 一 5 月 平均 归 一 化 植被 指数 空 
间 分 布 
Fig.5 
from 2005 to 2006 这 kuitufd river bas 盾 


Spatial distribution of averag@ NDVI in April and May 
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图 7 流域 土壤 含水 量 SRWC 值 空间 分 布 
Fig.7 Soil moisture spatial distribution of study area 
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图 8 土壤 含水 率 实 测 值 与 模拟 值 相关 分 析 
Fig.8 “Correlation analysis of soil moisture between measured 


and simulated value 
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(AS) 产 生地 表 径 流 为 边界 条 件 ,筛选 出 参数 率 定 期 112 组 、 验 证 期 58 组 符合 模型 条 件 的 降水 事件 。 再 利用 
聚 类 分 析 法 将 参数 率 定期 112 组 值 进行 聚 类 ,计算 出 30 组 聚 类 中 心 A 值 及 与 之 对 应 的 聚 类 成 员 ,将 这 30 
组 数据 作为 流域 初 损 率 (入 ) 取 值 库 ,参数 率 定期 径流 模拟 时 ,将 分 类 好 的 30 组 聚 类 成 员 分 别 带 人 对 应 的 聚 类 
中 心 A 值 进行 模拟 即 可 ;验证 期 径流 模拟 时 ,在 取 值 库 中 选择 与 其 计算 A 值 最 为 接近 的 聚 类 中 心 A 值 进行 模 
拟 即 可 。 对 干旱 区 大 尺度 流域 而 言 ,A 取 值 库 的 建立 较 传 统 单一 取 值 法 模拟 精度 更 高 。 


4 模拟 结果 检验 与 分 析 


本 文 利 用 改进 后 SCS 模型 对 奎屯 河流 域 参 数 率 定 期 (2005 一 2006 年 4.5 月 )112 场 降水 事件 和 验证 期 
(2007 年 4、5 月 )58 场 降水 事件 进行 径流 模拟 ,模拟 与 实测 径流 结果 如 图 9 所 示 。 
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9” 奢 屯 河 流域 参数 率 定期 改进 模型 模拟 值 与 实测 值 及 同期 降水 量 分 布 


Fig.9 Compared of observed and simulated runoff by calibration periods 


图 9、 图 10 分 别 绘制 了 模型 参数 率 定期 和 验证 期 模拟 值 与 实测 值 结果 。 从 模拟 结果 来 看 ,两 个 时 期 模型 
参数 算法 欧 进 后 的 模拟 值 与 实测 值 总 体 变化 趋势 一 致 , 且 与 流域 降水 量 分 布 呈正 相关 。 关 于 模拟 结果 的 精度 
检验 ,本 文选 取 汪 "Nash-Sutcliffe 效率 系数 和 相对 误差 系数 RE 来 反映 模型 模拟 值 和 实测 值 的 拟 合 程度 及 误差 
率 ( 表 2) 。Nash-Sttcliffe 效率 系数 经 常 被 用 来 作为 水 文 模型 的 效率 评价 指标 ,其 变化 范围 为 -ow 到 1。1 值 对 
应 于 实测 值 和 模拟 值 的 完美 匹配 ;效率 为 0 表明 该 模型 模拟 结果 等 同 于 实测 值 的 均值 系列 , 当 NSE 大 于 0.75 
时 本 认为 模型 模拟 结果 好 '] ,对 比分 析 模型 参数 率 定 期 和 验证 期 NSE 值 发 现 ,参数 率 定 期 模拟 值 与 实测 值 
匹配 程度 较 高 ,模拟 结果 效果 很 好 ;验证 期 NSE 值 略 低 ,说明 验 证 期 模拟 精度 还 有 竺 进一步 提高 。 相 对 误差 
系数 (RE) 是 最 能 反映 模拟 值 可 信 度 的 , 值 越 接 近 于 0, 则 说 明 模 拟 值 越 接近 实测 值 。 模 型 参数 率 定期 和 验证 
期 相对 误差 系数 RE 分 别 为 0.7% 和 -1.3% ,以 RE 低 于 +30% 为 标准 ; ,检验 模拟 结果 误差 范围 在 可 利用 
范围 。 

分 析 误 差 产 生 的 来 源 主 要 包括 以 下 几 个 方面 : 

(1) 计 算 $ 值 产生 的 误差 。 模 型 $ 值 的 决定 性 变量 是 土壤 含水 量 和 土壤 饱和 含水 量 。 土 壤 饱 和 含水 量 
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图 10 奎屯 河流 域 验证 期 改进 模型 模拟 值 与 实测 值 及 同期 降水 量 分 布 
Fig.10 Compared of observed and simulated runoff by validation periods 


取 最 大 值 1 ,会 导致 值 偏 大 引起 模拟 值 偏 小 的 误差 。 因 此 ,增加 土壤 舍 水 量 和 土壤 饱和 含水 量 野外 实测 数 
据 量 可 有 效 提高 5 值 的 精度 。 但 是 由 于 研究 流域 可 达 性 差 ,探讨 如 何 提 高 遥感 大 尺度 观测 数据 精度 将 会 成 为 
干旱 区 大 流域 径流 模拟 的 关键 技术 。 


表 2 参数 算法 改进 后 模型 模拟 值 效 率 系数 NSE、 相 关系 数 r 和 相对 误差 系数 RE 


Table2 NSE( Nash-Sutcliffe \efficiéncy coefficient) 、r ( correlation coefficient ) and RE ( relative error coefficient ) for simulation performance 


assessments 
评价 指标 参数 率 定期 验证 期 
Assessment indexes Calibration periods Validation periods 
纳什 效率 系数 NasisSutcliffe efficiency 0.92 0.64 
相对 误差 Relative errow/ % 0.7 -1.3 


(2) 计算 降水 量 造成 的 误差 。 基 于 流域 特殊 径流 补给 方式 ,改进 模型 中 的 降水 量 参数 被 修正 为 积 雪 消 融 
水 当量 和 实测 降雨 量 之 和 。 其 中 , 积 雪 消融 水 当量 的 决定 性 变量 是 温度 ,而 流域 内 唯一 的 观测 站 气温 数据 不 
居 以 代表 流域 平均 气温 ,而 积 雪 消融 水 当量 对 气温 变化 较为 敏感 ,流域 温度 每 变化 +1%C , 积 雪 消融 水 当量 变 
化 3.1 mm, 如 何 获 取 准 确 且 具有 代表 性 的 流域 气温 数据 也 是 提高 模拟 精度 的 关键 。 增 加 实测 观测 站 点 的 布 
设 , 受 到 地 理 位 置 和 经 费 等 问题 的 限制 ;就 目前 已 有 技术 而 言 ,可 利用 MODIS 日 数据 产品 提取 流域 温度 图 , 获 
取 流 域 平均 温度 ,从 而 实现 大 尺度 流域 气温 空间 数据 的 获取 。 

(3) 流 域 地 处 山区 ,海拔 高 ,山区 面积 大 ,分 布 有 冰川 ,流域 基础 地 理 数 据 和 野外 实测 数据 获取 十 分 困难 。 
因此 ,长 时 间 序 列 野 外 数据 的 缺乏 也 是 造成 模拟 误差 的 重要 方面 。 

此 ,在 今后 的 研究 中 ,探索 卫星 数据 和 遥感 技术 在 水 文 模型 中 的 应 用 ,增加 野外 实测 数据 量 等 都 是 提高 
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径流 模拟 精度 的 有 效 保障 。 


结论 


un 


(1) 为 了 使 SCS 模型 适用 于 奎屯 河流 域 径流 补给 特征 ,引入 度 -日 模型 修正 降水 量 参数 ,有 效 解决 了 因 降 
雨量 稀少 导致 模型 在 干旱 半 干 旱地 区 大 尺度 流域 无 法 应 用 的 瓶颈 。 

(2) 通 过 MODIS 卫星 数据 结合 7.-VI 特征 空间 理论 反 演 流域 土壤 含水 量 的 方法 ,集成 了 遥感 信息 空间 数 
据 的 时 空 分 布 特性 ,间接 实现 了 SCS 模型 参数 由 点 状 数据 向 面 状 数据 的 转换 。 

(3) 原 始 SCS 模型 计算 土壤 持 水 量 ($) 参 数 涉及 土地 利用 方式 等 5 个 参数 及 CN 值 对 照 表 ,利用 遥感 数 
据 计 算 $ 值 ,具有 涉及 参数 少 且 参数 易 获取 的 优势 ,可 有 效 解决 无 资料 地 区 模型 应 用 数据 缺乏 的 难题 。 

(4) 初 损 率 (A) 作 为 SCS 模型 重要 的 输入 变量 ,具有 时 空 变化 特性 ,采用 聚 类 分 析 法 建立 流域 入 取 值 库 
的 方法 ,使 模型 在 描述 流域 水 循环 过 程 时 更 接近 实际 ,其 模拟 结果 较 取 唯一 值 的 传统 方法 模拟 精度 有 较 大 幅 
提高 ,使 模拟 过 程 也 更 加 合理 。 

(5 ) 参 数 算法 改进 后 SCS 模型 参数 率 定 期 验证 期 NES(nash-sutcliffe efficieficy) 和 RE (relative error) 分 
别 为 0.92,0.64 和 0.7% ,-1.3% 。 参 数 改 进 算法 拓展 了 SCS 模型 的 适用 范围 ,可 为 其 他 干旱 半 干 旱地 区 下 热 
面条 件 相 似 流域 的 径流 模拟 提供 参考 。 
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